Mixed phases in U{N) superconductivity 
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We consider a general class of type-II superconductors which are described by TV complex order 
parameters (OP), with an overall U{N) symmetry for the energy functional. In the lowest energy 
state, each OP adopts its own particular lattice structure, making the spatial variation of the overall 
condensate very complex. We present detailed results for the case N = 2. We also treat the limit 
N oo, which allows us to resolve a recent controversy, by showing that all previous studies are 
fundamentally flawed. 



PACS: 64.60.Fr, 74.20.D 

I. INTRODUCTION 

It has been appreciated for many years that enlarg- 
ing the symmetry group of a given model can yield both 
new physical insights, and the possibility of exact calcu- 
lations. A familiar example is the generalization of the 
Heisenberg model of ferromagnetism to the 0{N) model 
[0 , which allows an exact analysis of the critical proper- 
ties for N ^ oo. The large- iV extension of the non- linear 
CT- model Q has been similarly fruitful Q. The idea has 
also been applied recently to non-equilibrium problems 
in turbulence phase ordering and interface rough- 
ening 1^. The implicit assumption involved is that the 
gross physical features are insensitive to the size of the 
symmetry group - an assumption which is largely borne 
out in the above applications. 

The first main point of this article is to emphasize that 
there is a wide class of models whose physics does not 
change smoothly on enlarging the symmetry group - in 
the language of critical phenomena, these are systems 
whose ordered phase is spatially inhomogeneous. We 
shall concentrate exclusively on one such system; namely, 
the type-II superconductor, which condenses into the so- 
called mixed-state - a periodic array of flux lines 0. 
However, non-trivial effects will be apparent in the other 
systems within this class (which may be fairly well char- 
acterized by having a wavelength selection mechanism in 
their ordered phase.) 

It is our intention to highlight the main physical 
changes which occur when one generalizes the Landau- 
Ginzburg theory of superconductors to a U{N) theory 
involving N complex order parameters (OP). The model 
is not new - in fact it has been studied many times in 
the past two decades or so [^-|T4[] - however, the physics 
we shall discuss here appears to have been overlooked in 
previous studies |10-13|. Although many of our conclu- 
sions apply to arbitrary N, such is the richness of these 
systems, we shall discuss in detail only the cases N = 2 
and N ~* oo. 

Our use of the U{N) symmetry is for convenience only, 
as it is the smoothest continuation of the model away 
from = 1. Clearly for higher values of N there are an 



increasingly large number of possible model Hamiltoni- 
ans which one may construct from the N OP's. Within 
the field of unconventional superconductivity (as evi- 
denced in the heavy- Fermion system UPta), the com- 
monly adopted model is an A^ = 2 Ginzburg-Landau the- 
ory with a reduced symmetry The richness of the 
phase diagram is well-documented, with the Abrikosov 
lattice giving way to exotic structures such as fractionally 
quantized vortices, and flux textures. An application for 
A^ = 3 can be found in rotating He^. Near to the normal- 
Al transition, the (18 component) OP may be reduced to 
three coupled complex fields Q , with the external angu- 
lar velocity playing the role of an applied magnetic field. 
These examples serve to show that the strategy of enlarg- 
ing the symmetry group of systems whose condensate is 
spatially inhomogeneous is a 'double-edged sword'. The 
advantage is that one finds an increasingly rich variety 
of phases as one increases A^. The disadvantage is that 
one thereby loses immediate physical contact with the 
original model of interest (A^ = !)• 

The remainder of this article is dedicated to exhibit- 
ing the strengths and weaknesses of the U{N) model, 
in terms of its relation to conventional superconductiv- 
ity. In section II, we define the model via the Landau- 
Ginzburg energy functional, and derive the mean field 
theory for arbitrary N. In section III we study in detail 
the mean-field solution for N = 2 and demonstrate that 
the Abrikosov state is unstable, giving way to a state of 
interlocked centered rectangular lattices. In section IV 
we review the large-A^ limit for this model. All previ- 
ous studies pO|-p^ have assumed an Abrikosov state for 
one condensed mode, and treated the remaining (N — 1) 
modes as massless. We show that this state is unstable, 
and that the true ordered state corresponds to a compli- 
cated structure in which each OP condenses into a pe- 
riodic state, such that the overall condensate density is 
spatially constant. This leads to an identification of the 
transition with that of the 0{2N) model of ferromag- 
netism (in the large-A" limit) in two fewer dimensions. 
An immediate consequence of this is that the transition is 
continuous and the lower critical dimension of the system 
is di = 4. We also discuss the subtleties of commuting 
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the limit N ^ oo with the thermodynamic limit. It is 
of note that the large- limit may be solved with no ap- 
proximations - we adopt neither the London limit (gauge 
field fluctuations only) nor the lowest Landau level (LLL) 
approximation (OP fluctuations only). We end with our 
conclusions in section V. 



II. FORMULATION OF THE MODEL AND 
MEAN FIELD THEORY 

We define the model via the Landau-Ginzburg energy 
functional which has the generic form (setting h — 
c=l) 

n[^,,A]= J d3a;{^[(l/2m*)|DV',|' + a|V'.|'] 

i 

+ (/3/2)^|V.n^,f + (l/2/.o)(H-Vx A)2} , (1) 

where D = — zV — e*A, {ipi} are a set of N complex order 
parameters, A is the vector potential, H is the external 
magnetic field, e* is the effective charge, and m*, a and /3 
are phenomenological constants. It is useful to introduce 
dimensionless units , such that the energy functional 
(in the ordered phase) takes the form 

i 

+ (l/2)^|V',nV',f + (H-VxA)2} , (2) 

where D — (l/iK)V — A. The parameter n > 1/^/2 is the 
ratio of the London penetration depth to the coherence 
length. In these units the critical external field Hc2 — k. 
When generalizing the above expressions to d dimensions, 
one should imagine the external field directed along {d — 
2) dimensions, such that it is transverse to the (x, y) 
plane in which an Abrikosov-like lattice structure may 
form. 

The simplest analysis one can make of this system is 
mean field theory which amounts to approximating the 
free energy by the saddle point value of Ti. This may 
be obtained by solving explicitly for the classical fields 
which are solutions of the Landau-Ginzburg equations 

DV. = (i-EI^^-|')'/'^ (3) 

j 

V X b = (l/2)(v*DV'-h^/'D*V'*) , (4) 

where b = H — VxAis the microscopic magnetic field. 

Although it is not known how to solve these equations 
in general, an exact analysis is possible very close to the 
critical field Hc2- We shall briefly outline the solution in 
the case of iV = 1, and then generalize the solution to 



arbitrary N. Further details may be found in refs. [f7|,pT|. 
Near to Hc2 the OP is small, so that the first equation 
may be linearized, and at this order the gauge field is 
just given by Aq = Hc2{0, x,0±). Thus, the OP may 
be expressed in terms of Landau levels [^. In fact, the 
value of Hc2 itself is determined by associating criticality 
with the LLL eigenvalue. For our purposes, the key point 
is that the first equation may be rearranged in the form 



where 



D+D.^P = , 



(5) 



(6) 



and we have used the fact that the OP components are 
constant in the (d—2) other directions. Thus the ground 
state 'wavef unctions' are characterized by the identity 



D_ijj = 



(7) 



These functions are the LLL, and have a degeneracy pro- 
portional to the area of the system (transverse to the 
applied magnetic field). One is free to construct different 
sets of functions from the LLL's. The most elegant is due 
to Eilenberger |]l9| . The Eilenberger basis functions span 
the space of the LLL, yet each basis function is a dou- 
bly periodic function (essentially a particular Abrikosov 
lattice). 

Given the above properties of the OP near the upper 
critical field, it is possible to manipulate the second mean 
field equation such that the r.h.s. is the curl of a vector. 
In this way, one may integrate the equation to obtain 
Abrikosov's first fundamental identity 



b{r) = H - {l/2K)\i;{r)\' 



(8) 



We see that the magnetic field is reduced within the sam- 
ple by an amount proportional to the condensate density. 

The scale of the condensate is set by Abrikosov's sec- 
ond identity. This is a little more difficult to prove , 
but is essentially obtained by spatially averaging the first 
equation and calculating the first non-trivial corrections 
to eq.(|7|). The result is 

- Hm^) + ((1/2^2) _ ^)^|^^|4^ ^ ^ 

where (...) denotes a spatial average. 

Substituting this last relation into eq. (|^) and spatially 
averaging, we find 



B^(b) =H 



(k-H) 



(2«2 _ i)/3^ ' 
where the Abrikosov ratio is given by 

(3A = {m/m'r>i. 



(10) 



(11) 



One can translate this result in terms of the free energy, 
which one finds to be proportional to —I/Pa, indicating 
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that the system condenses into a state which minimizes 
Pa under the constraint that it is an Eilenberger func- 
tion. This state turns out to be the well-known triangular 
lattice, for which (3a — 1.1596. 

Now we come to the general U{N) case as described 
by eqs.(||) and (|^). It is easy to show that most of the 
previous analysis follows through in a trivial way. Each 
OP component must satisfy D-ipi — (although it is 
important to remember that a given component may be 
zero). The first Abrikosov identity is generalized to 



b{r) = H - {l/2K)J2\Mr)\' 



There also exist N relations setting the relative scales of 
the OP components. These have the form 

(l/^)(^ - H)m') + ((1/2^2) _ 1) ^(|V>.nv.,f ) = . 

j 

(13) 

This allows us to spatially average eq.(p^) to find the 
relation between the magnetic flux density and the mag- 
netic field: 



B = H 



{n-H) 



(2^2 _ i)(3^[N) 



where the generalized Abrikosov ratio is given by 



(E.dV'.P))' 



> 1 



(15) 



Again, one may show that the minimal free energy is 
to be found by minimizing /3g, with the constraint that 
each condensed OP component is either zero, or an Eilen- 
berger function. For general N one soon appreciates that 
such a minimization is non-trivial as the N components 
jostle for favorable positions and lattice structures within 
the 'primitive cell'. For this reason we concentrate in the 
next section on the simplest case. 



III. ANALYSIS OF THE CASE N = 2 

We are now only concerned with two OP's ipi and ■02- 
By reducing the U{2) symmetry we could make contact 
with the two-OP models of UPt3 [|5| for which a similar 
mean field analysis has been performed pof . Our purpose 
here is to exemplify the physics of these multi-component 
systems - this will be of great benefit in our analysis of 
the large- limit in the following section. 

A few words concerning the Eilenberger basis are 
required at this point. As mentioned before, the Eilen- 
berger functions (/)(r|ro) satisfy Z3_0 = 0. The amplitude 
of is a doubly periodic function, with a fundamental cell 
scaled so as to have unit length in the x direction, and a 



periodicity vector (C, rj), where rj is fixed by the condition 
of flux quantization. The label Tq simply fixes the lattice 
position in space. The functions are normalized such that 
{\4>\^) — 1. The functions also have the symmetry prop- 
erty 0(r|ro) = exp[27ri(yo/??)a:]0(r + ro|0). This allows 
one to recast the normalization condition as a complete- 
ness relation; 



dVo|(/.(r|ro)|2 = 77 . 



(16) 



(12) It is useful to define the integrals 



/(ri,r2) 



dV|0(r|ri)n0(r|r2) 



(17) 



Returning to our expression for j3g in the case N 
we have to minimize the expression 



/3s (2) = 



(|V2|^)+2(|V'l|2|02 



((1^1 



(IV'2P))^ 



(18) 



One may show that a minimum may only exist for ei- 
ther one of the OP's being zero, or else for each OP 
to be equivalent (up to a relative spatial shift). In 
(14) the former case the value of fig will be the Abrikosov 
value (~ 1.1596). To investigate the latter case, we set 
V"! = A(/>(r|0) and V'2 = A(f){v\vQ) which reduces the task 
of minimizing (3g to that of finding the lattice type and 
the relative shift ro which minimize 



/?g(2) = 7(0,0) + 7(0, ro) 



(19) 



We have numerically investigated the above expres- 
sion, restricting our search to lattices within the class 
of centered rectangular structures (which corresponds to 
choosing C, — 1/2). This class includes square and trian- 
gular lattices. Somewhat surprisingly the minimal energy 
solution corresponds to each OP component adopting a 
lattice with a primitive cell with opening angle 9 equal 
to 15°. Figures 1 and 2 show contour plots of the two 
OP components. Regions of lighter shade correspond to 
higher values of the OP. It is more illuminating to plot 
the overall condensate p -I- |(/)2p, as is shown in Fig. 3. 
We then see a surprisingly rich structure. Maxima of the 
condensate correspond to regions of minimal magnetic 
flux and vice versa. The energy for this arrangement 
corresponds to the value /3g(2) ~ 1.0062, which is very 
close to the lower bound of unity. 

Our initial guess was that the OP components would 
adopt square lattices {0 = 45°), as this arrangement has 
a higher symmetry. We plot the components and the 
overall condensate in Figs. 4, 5 and 6 for this case. The 
energy of the square lattice arrangement is only fraction- 
ally higher with [ig{2) ~ 1.0075. 



3 



FIG. 1. 101 1^ for minimum energy lattice. FIG. 4. |0i|^ for square lattice. 




FIG. 3. Condensate + |<^2|' for minimum energy lat- p^Q q Condensate \4>i\'' + |02|' for square lattice. 



In Fig. 7 we show the generalized Abrikosov number 
as a function of angle. Interestingly, the triangular lat- 
tices correspond to maximal values of (3g. In principle, 
there may be an even lower energy OP arrangement cor- 
responding to a lattice with an oblique primitive cell lying 
outside the class of centered rectangular structures. This 
is a large space within which to search for minima, and 
we have not pursued this possibility. 
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FIG. 7. The function /3g(2) plotted against 6. The function 
is symmetric about 45°. 

From this analysis, we see that in general the U{N) 
models will condense into complicated structures in 
which each component adopts a particular lattice con- 
figuration, such that the overall condensate density is as 
smoothly varying as possible. The calculation of these 
structures becomes increasingly difficult as one increases 
N, and we shall desist from any explicit calculations In 
the next section we turn to the most important part of 
the article - namely the explicit solution of the U{N) 
model in the large- limit. 

IV. THE LARGE-Af LIMIT 

As hinted at in the Introduction, one of the more com- 
pelling reasons for generalizing models to higher sym- 
metry groups is to allow exact solutions in the large- A^ 
limit. This has proven to be a very useful tool in many 
contexts . In the present case, the large- A^ limit of 

the U{N) model of type-II superconductors (in an exter- 
nal field) has been examined by several authors [^-^. 
[One should draw a clear line between these calculations, 
and those concentrating on the zero field case which 
is of interest in liquid crystals and also as an exemplifica- 
tion of a fiuctuation induced first-order phase transition 
Originally studied in 1985 [|o|, it was found that 
below six dimensions, the mean-field continuous transi- 
tion of Abrikosov becomes first-order - all the way down 
to a lower critical dimension of d/ = 2. This calculation 



relied upon a proof by contradiction, as an explicit solu- 
tion of the model is extremely difficult. An independent 
study was made in 1995 using an Ansatz to solve 
the model. It was found that the transition was contin- 
uous below six dimensions, and that di — 4. This result 
was challenged on the grounds that the Ansatz was 
physically unreasonable and that the original calculation 
pO{ was definitely correct. A spirited response |]l3| was 
made, defending the Ansatz, but admitting that there 
was no obvious fiaw in the original proof by contradiction 
pO| . This is quite an extraordinary situation regarding a 
model whose raison d'etre is its solvability! 

We shall resolve this state of confusion using a very 
simple physical idea that has emerged from the previ- 
ous two sections. Each of the previous studies follows 
the conventional route of integrating out (A^ — 1) of the 
OP components, regarding them as Goldstone modes. 
This is precisely what one does in the conventional 0{N) 
model of ferromagnetism for instance ||l|,^ . The reason 
it is done there is that the condensed OP is spatially ho- 
mogeneous, and it makes good sense to globally rotate 
the OP's such that the condensate exists exclusively in 
one component (denoted as longitudinal), and to treat 
the transverse components as a source of massless fiuc- 
tuations. Such a course of action is inadmissible in the 
present model of U{N) superconductivity. Each compo- 
nent of the condensate is spatially inhomogeneous and 
no global rotation is possible to make (A^ — 1) Goldstone 
modes. One has to treat all the modes as potentially 
condensed. For finite A^ this would lead to a totally in- 
tractable problem. However, for A^ — s- oo we can take ad- 
vantage of the completeness relation of Eilenberger func- 
tions to demonstrate that although each OP component 
is a doubly periodic function, the overall condensate is 
a constant. This drastically simplifies the analysis, and 
one may show that the solution is self-consistent when 
one takes into account both higher Landau levels, and 
gauge field fluctuations. 

Before turning to the self-consistent large- A^ limit, we 
shall obtain the main idea by studying the limit of large 
A^ at mean field level. The task is to minimize (3g{N) as 
given in eg. ([l5|) . We see that by choosing each compo- 
nent to lie in a different Eilenberger state, the sums over 
components tend (in the limit of large A^) to integrals 
over the basis label Tq. We may then invoke the com- 
pleteness relation given in eq.(|l6|) which reduces these 
component sums to constants. Hence the spatial aver- 
ages are trivial, and we see that /3g is saturated at its 
lower bound of unity. 

It is a remarkable fact that this particular solution 
actually solves the mean field equations throughout the 
mixed phase, and not just in the vicinity of the 11^2 line. 
One can see this straightforwardly: since J^i iV'iP is a 
constant, the first mean-field equation is exactly solved 
by the Eilenberger functions. The higher Landau levels 
always remain hard modes as one decreases the temper- 
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ature, so never contribute to the condensed OP. At low 
enough temperatures, the Meissner transition will occur, 
and the condensate will saturate at the value 

2kH. 

Two interesting points follow. In the previous section 
we found that even for two OP components, the con- 
densed state was very smooth (in terms of the overall 
condensate). In that sense it is already very similar to 
the large- limit in which the condensate is exactly con- 
stant. Thus we expect very similar quantitative physics 
as one increases N in the range [2,oo]. Note also the 
huge degeneracy of lattice structures underlying this so- 
lution. Although each OP component is restricted to be 
an Eilenberger function with a particular shape of unit 
cell, there is no energetic selection of that unit cell for 
N oo, since the completeness relation is independent 
of this. Presumably the ground-state OP configurations 
for large but finite N provide a means of determining the 
large- A^ configurations, by smoothly continuing A^ — > cxd. 

The self-consistent large- A^ treatment is a non-trivial 
extension of the above calculation, as it takes into ac- 
count fluctuations about mean-field theory (albeit in a 
rather crude fashion.) We shall find that the main char- 
acteristics of the mean field state are stable to these fluc- 
tuations for d > di = A. For 4 < d < 6 the fluctuations 
give rise to cZ-dependent exponents, which cross over to 
mean-field values above six dimensions. In fact, from the 
structure of the theory, one sees that the results are iden- 
tical to those of the large-A^ limit of the 0{2N) model 
of ferromagnetism in two dimensions fewer. We can un- 
derstand this as the correlations of the OP's are frozen 
in the {x, y) plane due to each OP component having 
formed a lattice state, so only transverse fluctuations (in 
the remaining [d — 2) dimensions) can become critical. 

[We note here that these results are similar to those of 
ref. but we must stress that the physical condensed 
state is completely different in the two cases. In ref. 
there is no explicit transverse scale as only one OP com- 
ponent is condensed, and is taken ad hoc to be constant 
(which is not even a solution of the saddle-point equa- 
tions ||l^ ); whereas in our solution, each OP component 
has condensed into a lattice state with its own magnetic 
length scale built in. Real thermal fluctuations about 
these two states will be of totally different natures.] 

The large-A^ limit is often derived diagrammatically 
[0,^. When one can describe the physics in terms of 
a longitudinal mode and (A^ — 1) transverse modes, this 
approach is particularly transparent. However, in the 
present case we have N condensed modes, so it is bet- 
ter to use an alternative method; namely to introduce 
an auxiliary field x (via a Hubbard-Stratonovich trans- 
formation) which will allow us to make the large- A^ limit 
explicit. 

In the presence of fluctuations, we must take a step 
back from eq. (^) and consider the partition function Z ~ 
J VAVip, exp[-7^'], where 



W =H- I d'^r^(J*V^ 



+ c.c) 



(20) 



The source terms Ji are added in order to derive the 
equation of state in the ordered phase. It will turn out 
that each Ji is proportional to an Eilenberger function. 

Introducing the field x allows us to rewrite the parti- 
tion function in the form 



Z[J] 



VAVxVip^ exp[-7^"] 



(21) 



where 

n"= j d'^r{Nx^/W + N{n - V X A)2 

i 

We have scaled the magnetic field and the vector poten- 
tial so as to extract a clean factor of A^ in the first two 
terms [|l^. (This entails the rescalings e* — > e* / ^JN , and 
(3 P/N for consistency). 

As we have already indicated, the self-consistent large- 
A^ limit for this problem constitutes a formidable analytic 
challenge. In fact, it is not possible to solve the problem 
without resort to some external resource, whether it be 
an Ansatz, or a piece of physical insight. We shall utilize 
the latter, thanks to the lessons we have learned both in 
the A = 2 case, and also in the mean field analysis of 
the large-A^ limit. Just to reiterate, at mean field level 
the OP components each condense into an Eilenberger 
state such that the overall condensate is spatially homo- 
geneous, and consequently the magnetic flux is spatially 
homogeneous also. To proceed we take the simplest pos- 
sible line. Namely, that the self-consistent treatment re- 
tains these features of spatial homogeneity, but that the 
fluctuations renormalize the mass a, resulting in a shift of 
Tc- Our task is to show that this is a consistent solution 
of the problem. The alternative Ansatz is to condense 
only one OP component [^o|-^3|. Rewriting the energy 
functional in terms of Eilenberger functions allows one to 
prove that such a state is energetically unstable |2^ . 

To select the physically motivated condensed state we 
must choose the source terms to force each component 
into a (spatially shifted) Eilenberger function. Thus we 
write Ji = u(j)i — u(j)(r\ri), where u is complex. The 
homogeneity of the magnetic field allows us to write B = 
5(0, 0, fj,). Also it is convenient to set t = r -|- ix- The 
saddle point value of x is purely imaginary, such that the 
effective mass t is purely real. The energy function now 
takes the explicit form 



H" = J d'^r{~N{t - rf/2f3 + N{H - Bf 

+ IIW'^^I' + (l/2m*)|DV'»P - u4>*^, ~ u<j),4,*] } , (2 
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where now A = B{0, x, 0±). Each individual OP compo- 
nent is decoupled, and may be associated with a partition 
function z{t,u) — J V^jji exp{~ f) where 



d'^r ['^*M^P^ - u(j,*4,, - udp,4,*] , (24) 



where M — —{l/2m*)TP + i, which has eigenvalues 

A.,. = ^+^+(-+2)^, (25) 

where the momentum k exists in the {d — 2) dimensions 
transverse to the (x,y) plane, and n e [0,oo] labels the 
Landau levels. 

We now write the condensed part of the OP explicitly, 
along with its fluctuation: -0^ = uxpi + tpi. The prefactor 
w is chosen to ensure that the fluctuation piece ipi has 
zero mean. Substituting this into eq. ( p^ ) and performing 
the integrals over the Eilcnbcrgcr functions, one finds 
that the terms linear in ij^i vanish so long as one chooses 
w = u/{t + e* B /2m*). In this case the energy functional 
for each component becomes 



V\u\ 



2{t + e*B/2m* 



+ / d'^r ^p*M^, , (26) 



where V is the volume of the system. 

Returning now to the energy functional, we have from 
eqs.(|2^) and (p6|), along with assumed spatial constancy 
of X and B, 



n" = - NV 



2/3 



(H-B) 



2{t + e*B/2m*) 



+ / d\ ^*Mi^, 



(27) 



The integrals over the fluctuation fields {"0^} are easily 
done, and we may re-exponentiate the resulting determi- 
nant to give the final result 



NV 



2(3 



-{H-Bf 



2{t + e*B/2m*) 



(Dl/A) 



log(Ak,„) 



(28) 



where Dl is the Landau degeneracy of each level, which 
is equal to e*BA/2Tr, where A is the area of the system 
in the (a;, y) plane. 

Now that we have a clean factor of throughout, we 
may use steepest descents to determine the self-consistent 
values of the auxiliary variable t , and also the magnetic 
fleld B. Differentiating H" with respect to t yields 



2{t 



e*B/2m*)2 

e*B r d'^-^k 
2n J (27r)'^-2 



(29) 



It is convenient to define = t + e* B /2m* , since ^ can 
be seen to play the role of the correlation length of the 
system. The relation between w and u then assumes the 
form w = u^^. One may then rewrite eq.(p9|) as 

^-^= r + e*B/2m* + (3\w\^/2 



e*B r d'^-^k ^ 



[fc2/2TO* + ne*B/m* + ^-2] 



(30) 



This equation encapsulates most of the information 
about the phase transition. Above the transition, we can 
set the condensate 'amplitude' w to zero in (^. The 
resulting equation defines the renormalized critical tem- 
perature Tc through the defining condition of criticality 
^ — + cxD. At this point the bare quantity r = T — is 
equal to Tc — T°. We therefore have the explicit shift as 



Tc = TO - e*B/2m 
e*B r d'^-'^k 



27r 



(27r) 



d-2 5Z 



[k'^/2m* + ne*B/m* 



(31) 



As usual, the fluctuations drive the critical temperature 
to a lower value. The Tc shift diverges for d < 4, sug- 
gesting the identiflcation of the lower critical dimension 
as di — 4. 

On the low temperature side of the transition we can 
remove the source field u, which leaves a non-zero con- 
densate w only if the correlation length is infinite. We 
therefore have from (^0|), an equation for the condensate 
amplitude: 



^r + e*B/2m*+f3\wf/2 



e*B 
~2^ 



d'^-^k 



(2^) 



d-2 ^ 



[fc2/2m* -I- ne*B/7 



(32) 



Comparing eqs.(|3l|) and ( p2D we find the exact relation 
|w|2 ^ 2{Tc-T)/l3 which immediately yields the OP 
exponent [3 = 1/2, and self-consistently confirms the ex- 
istence of a continuous transition. 

We can also identify the correlation length exponent 
by examining ( ^0|) as T \ Tc. Eliminating the bare 
critical temperature from ( |30|) using ( |3l]) , and evaluating 
the resulting integral for large ^ leads to the expression 



l3e*Bc2£.^-'^ = T -Tc 



(33) 



The constant ci has a contribution from all Landau levels 
bar the lowest. The fluctuation dominated term ~ ^^-d 
and arises solely from the LLL. As the correlation length 
diverges for T \ Tc, we see that the first term on the 
left hand side dominates for d > 6, whereas the second 
term dominates for 4 < rf < 6. This leads to the result 
^ ~ (T - Tc)-'', with 1/ = 1/2 for d > 6 (the mean field 
result) , and v — 1/ (d — 4) for 4 < d < 6 (confirming 
di = 4). 
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These results are identical to those obtained for the 
0{2N) model of ferromagnetism, but in two fewer di- 
mensions. This may be understood from examination of 
the self-consistent relation for ^ given in ( |30| ) . Apart from 
the sum over Landau levels, this equation is exactly that 
which would be obtained for a (d— 2)-dimensional 0{2N) 
model. The critical modes only exist in the {d — 2) di- 
mensions transverse to the Landau levels. The modes in 
the {x, y) plane contain the frozen length scale associated 
with the formation of the underlying lattice structure of 
the OP's. 

In the above expressions, we have left the value of the 
magnetic field B undetermined. However, this is given 
self-consistently by minimizing the energy functional in 
(^) with respect to B. We shall not write the expression 
explicitly, but it is noteworthy that the integral appear- 
ing from the fluctuations is strongly divergent and must 
be regularized by introducing some microscopic cut-off 
procedure (for d > 4), such as adding higher derivative 
terms not present in our original Landau-Ginzburg en- 
ergy functional. 

As a final remark in this section, we should point out 
that these results are insensitive to the order in which one 
takes the thermodynamic limit, and the limit iV — > oo. 
The above analysis has implicitly assumed the thermody- 
namic limit, which is the correct 'physical' choice. How- 
ever, had one taken the thermodynamic limit second, by 
fixing the number of vortices and then taking N oo, 
the calculation would have proceeded as before with one 
difference. The fluctuations in this case would origi- 
nate overwhelmingly from the transverse modes, since 
the number of distinct longitudinal modes is limited to 
the number of vortices. However the completeness rela- 
tion still holds, and the condensed state may be taken as 
spatially homogeneous, thereby enabling a calculation in 
the same spirit as that above. 

V. CONCLUSIONS 

The main point of this article is that special care must 
be taken when expanding the symmetry group of systems 
with spatially varying structures. We have concentrated 
on one class of such systems, namely type-H supercon- 
ductors in an external magnetic field. As mentioned in 
the Introduction, there are applications of such models 
to heavy- Fermion superconductors, and also to rotating 
superfluid Hea. We have found a number of interesting 
results connected with A^-component superconductors, 
whose free energy functional maintains a U{N) symme- 
try. In general we have seen that these systems adopt 
low-temperature configurations in which many OP com- 
ponents contribute by condensing into periodic structures 
- leading to a very rich structure for the overall conden- 
sate (and hence the magnetic fiux). 



We have examined the mean-field theory for the case 
A^ = 2 in some detail. The two OP components were 
found to condense into centered rectangular structures, 
with an opening angle of 15°. The two structures are 
shifted relative to one another in such a way that the 
overall condensate + |(/>2p has a surprisingly rich 

structure, as shown in fig. 3. The generalized Abrikosov 
ratio for this configuration is /3g(2) ~ 1.0062, almost 
saturating the lower bound of unity. Although systems 
with higher values of A'^ will adopt ever more complex 
structures, we were able to show in section IV that the 
mean-field theory in the limit N oo has a simplifying 
feature, due to the completeness relation of the periodic 
Eilenberger functions. Each OP component adopts an 
Eilenberger function, but the overall condensate has no 
traces of the periodicity, and is spatially a constant. This 
structure saturates the lower bound of l3g{oo). 

In the remainder of section IV we examined the large- 
A^ limit in more detail, by considering a treatment which 
includes fiuctuations self-consistently. This calculation 
has been attempted several times in the past [p^-|l3| , but 
the previous authors have always assumed that the sys- 
tem allows only one OP component to condense, which 
we have seen is generically false. Our main finding is that 
for d > 4, fluctuations do not disturb the main charac- 
teristics found in mean-field theory - namely a continu- 
ous transition into a spatially homogeneous condensate, 
composed of infinitely many OP components having con- 
densed into Eilenberger functions. [For d < 4, the mixed 
phase is destroyed entirely] The system is found to have 
the critical properties of the 0{2N) model of ferromag- 
netism p],p2|, but in two fewer dimensions. Thus, expo- 
nents maintain their mean-field values above d = 6, but 
take the values (3=1/2 and v = l/(d- 4) for 4 < d < 6. 
The solution we have found does not rely upon making 
the LLL approximation, or upon neglecting gauge field 
fluctuations. 

In this paper we have found that within mean-field 
theory, and also in the large- A^ limit, the U{N) model 
undergoes a continuous transition from the normal to 
the mixed phase. This is in contradiction to the results 
of some past works |^,|l^,^ , the latter two of which con- 
tain errors of principle. However, one may find a prece- 
dent for a continuous transition in our recent study [|l4| , 
in which a functional renormalization group (FRG) ap- 
proach was applied to the U{N) model via an expansion 
in e = 6 — d. In fact, the FRG study also predicted a 
mapping from the U{N) model to the 0{2N) model of 
ferromagnetism in two fewer dimensions, for N > 2. It 
would be interesting to probe this relationship further 
by extending the present self-consistent analysis of the 
U{N) model to finite A^. One of the more sophisticated 
means of achieving this would be via the use of the par- 
quet approximation p^ ], which includes corrections far 
beyond those of 0{l/N). 

Finally we would like to draw the reader's attention to 
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the fact that in the exactly solvable large- iV limit, the 
mechanism for the transition is the growing (phase) co- 
herence in the direction of the applied field, as the tem- 
perature is lowered. A theory of the "melting" transi- 
tion seen in high-temperature superconductors {N =1) 
has been recently given by one of us based on the 
idea that the apparent melting is just a consequence of 
crossover effects, when this phase coherence length scale 
(which is very rapidly growing in three dimensions), be- 
comes comparable to the dimensions of the system. 

TJN acknowledges financial support from the Engineer- 
ing and Physical Sciences Research Council. 
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